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Abstract 

We derive predictions for the Nadir angle (On) dependence of the ratio N^~ /N^+ of 
, the rates of the fi— and //"'" multi-GeV events, and for the fi~—fi~^ event rate asymmetry, 

\^ I ^f_i~fj.+ = [A(/x~) — A^(/i+)]/[A(/x~) + A(;U+)], in iron-magnetized calorimeter detectors 

0^ ■ (MINOS, INO, etc.) in the case of 3- neutrino oscillations of the atmospheric v^ and 

^^ , Uf^, driven by one neutrino mass squared difference, |Arn-|]^| ~ (2.0 — 3.0) x 10~^ eV^ 3> 

(^ I Attt-Ii- The asymmetry ^ + (the ratio AL-/A^+) is shown to be particularly sen- 

's!" ■ sitive to the Earth matter effects in the atmospheric neutrino oscillations, and thus 

^^ , to the values of sin^ ^13 and sin^ 623 , ^13 and 623 being the neutrino mixing angle 

limited by the CHOOZ and Palo Verde experiments and that responsible for the dom- 
inant atmospheric v^ — > Vr (i^/^ -^ ^r) oscillations. It is also very sensitive to the 
^H. type of neutrino mass spectrum which can be with normal (Aml^ > 0) or with in- 

r-| ! verted (Am^^ < 0) hierarchy. We find that for sin^ 6*23 ^ 0.50, sin^ 2^13 ^ 0.06 and 

|Ar?T,3^| = (2 — 3) X 10~^ eV'^, the Earth matter effects produce a relative difference bet- 
ween the integrated asymmetries ^^-^+ and ^^-„+ in the mantle (cos0„ = 0.30 — 0.84) 
and core (cos 9n = 0.84—1.0) bins, which is bigger in absolute value than approximately 
C^ I ~ 15%, can reach the values of (30 — 35)%, and thus can be sufficiently large to be ob- 

servable. The sign of the indicated asymmetry difference is anticorrelated with the sign 
of Arn^i. An observation of the Earth matter effects in the Nadir angle distribution 
of the asymmetry A^-^+ (ratio N^~ /N^+) would clearly indicate that sin^ 2^13 ^ 0.06 
and sin^ 623 ^ 0.50, and would lead to the determination of the sign of Hwn^i. 
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1 Introduction 

There has been a remarkable progress in the studies of neutrino oscillations in the last 
several years. The experiments with solar, atmospheric and reactor neutrinos ^1 12 El IH 
El El Ej have provided compelling evidences for the existence of neutrino oscillations driven 
by nonzero neutrino masses and neutrino mixing. Evidences for oscillations of neutrinos 
were obtained also in the first long baseline accelerator neutrino experiment K2K [S] . It was 
predicted already in 1967 9j that the existence of solar neutrino oscillations would cause a 
deficit of solar neutrinos detected on Earth. The hypothesis of solar neutrino oscillations, 
which in one variety or another were considered starting from the late 60's as the most natural 
explanation of the observed [H EI solar neutrino deficit (see, e.g., refs. 01111111111121), has 
received a convincing confirmation from the measurement of the solar neutrino flux through 
the neutral current reaction on deuterium by the SNO experiment jH |3| , and by the first 
results of the KamLAND experiment 0- The analysis of the solar neutrino data obtained 
by Homestake, SAGE, GALLEX/GNO, Super-Kamiokande and SNO experiments showed 
that the data favor the Large Mixing Angle (LMA) MSW solution of the solar neutrino 
problem (see, e.g., ref. Q). The first results of the KamLAND reactor experiment jZj have 
confirmed (under the very plausible assumption of CPT-invariance) the LMA MSW solution, 
establishing it essentially as a unique solution of the solar neutrino problem. 

The latest addition to this magnificent effort is the evidence presented recently by the 
Super-Kamiokande (SK) collaboration for an "oscillation dip" in the L/ii^— dependence, of 
the (essentially multi-GeV) /i— like atmospheric neutrino events ^ jE], L and E being the 
distance traveled by neutrinos and the neutrino energy. As is well known, the SK atmospheric 
neutrino data is best described in terms of dominant 2-neutrino v^ —>■ v^ [u^ — »• z/,-) vacuum 
oscillations with maximal mixing, sin^ 2^23 — 1- The observed dip is predicted due to the 
oscillatory dependence of the Ufj, -^ Ur and u^ -^ Ur oscillation probabilities, ^(t'/i — ^ i^t) — 
Pi^fi -^ ^t), on L/E. The dip in the observed L/E distribution corresponds to the first 
oscillation minimum of the z/^ (z/^) survival probability, P{i^^ — ^ z/^) = 1 — P{i^fi — > i^r), 
as L/E increases starting from values for which Am'^^L/{2E) <^ 1 and -P(z^^ — > z/^) = 1, 
Amg^ being the neutrino mass squared difference responsible for the atmospheric z/^ and u^ 
oscillations. This beautiful result represents the first ever observation of a direct effect of 
the oscillatory dependence on L/E of the probability of neutrino oscillations in vacuum. 

The interpretation of the solar and atmospheric neutrino, and of KamLAND data in terms 
of neutrino oscillations requires the existence of 3- neutrino mixing in the weak charged lepton 
current (see, e.g., ref. [HI): 

3 

r^iL = Y.UiJ^J^ ■ (1) 

i=i 

Here un, I = e,fi,T, are the three left-handed flavor neutrino fields, z/jl is the left-handed 
field of the neutrino Uj having a mass rrij and U is the Pontecorvo-Maki-Nakagawa-Sakata 



"'^The sample used in the analysis of the L/E dependence consists of /i— like events for which the relative 
uncertainty in the experimental determination of the L/E ratio does not exceed 70%. 



(PMNS) neutrino mixing matrix [T3] . 
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( C12C13 S12C13 

-S12C23 - Ci2S23Si3e*'' C12C23 - Si2S23Si3e*'' 
V '^12'S23 "" Ci2C23Si3e* — C12S23 — Si2C23Si3e* 




where we have used a standard parametrization of U with the usual notations, Sij = sin Oij 
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jj = cos6ij, and 6 is the Dirac CP-violation phase ^. If one identifies Am2i > and A?7i|^ 
with the neutrino mass squared differences which drive the solar and atmospheric neutrino 
oscillations, the data suggest that lATTig^^l ^ Att^Ii- In this case 612 and 623, represent the 
neutrino mixing angles responsible for the solar and the dominant atmospheric neutrino 
oscillations, 612, ^23, while ^13 is the angle limited by the data from the CHOOZ and Palo 
Verde experiments [121 120] • 

The 3-neutrino oscillations of the solar z/g depend in the case of interest, | Amg^^l ^ ArTi^i, 
not only on Am2i and 612, but also on 613. A combined 3-neutrino oscillation analysis of 
the solar neutrino, CHOOZ and KamLAND data showed |2I] that for sin^6'i3 ^ 0.05 the 
allowed ranges of the solar neutrino oscillation parameters do not differ substantially from 
those derived in the 2-neutrino oscillation analyzes (see, e.g., refs. [3 123)- A description of 
the indicated data in terms of z/g — > z/^.^ and z/g -^ z/^^ oscillations is possible (at 99.73% 
C.L.) for sin2^i3 < 0.075. The data favor the LMA-I MSW solution (see, e.g., refs. |2III22I) 
with Am|i = 7.2 x 10~^ eV^ and sin2 6'i2 = 0.30. The LMA-II solution, corresponding to 
Am^i = 1.5 X 10~^ eV^ and approximately the same value of sin^ 612, is severely constrained 
by the fact [SB] that the ratio of the rates of the CC and NC reactions on deuterium, measured 
with a relatively high precision in SNO during the salt phase of the experiment [Sj, turned 
out to be definitely smaller than 0.50. This solution is currently allowed by the data only at 
99.13% C.L. |21i. 

The preliminary results of the most recent improved analysis of the SK atmospheric 
neutrino data, performed by the SK collaboration, gave at 90% C.L. [21] 

1.3 X 10"^ eV^ < \Aml^\ < 3.1 x 10"^ eV^ , 0.90 < sin^ 2^23 < 1-0, (3) 

with best fit values |Am^^| = 2.0 x 10"^ eV^ and sin^ 26'a = 1-0. Adding the K2K data [H], 
the authors [2H] find the same results for sin^2^A5 the same lArrig^l best fit value and 

1.55 X 10"^ eV^ < lAmiil < 2.60 x 10"^ eV^ , 90% C.L. (4) 

Earlier combined analyzes of the SK atmospheric neutrino and the K2K data produced 
somewhat larger values of |Am|^|: the best fit value and the 90% C.L. allowed interval 
of values, found, e.g., in ref. [26| read lArrig^l = 2.6 x 10~^ eV and 2.0 x 10"^ eV^ ^ 
lAmg^l ^ 3.2 X 10~^ eV^. Finally, the values of jAmg^l and sin^26'23, deduced from the SK 
analysis of the L/E dependence of the observed /x— like atmospheric neutrino events [T3] . 
are comfortably compatible with the values obtained in the other analyzes: |Am|]^| = 2.4 x 
10-3 eV, sin2 2^23 = 1 (best fit), and 1.9 x 10"=^ eV^ < |Am|i| < 3.0 x 10"=^ eV^ 0.90 < 
sin^ 26*23 < 1-0 (90% C.L.). We will use in our further discussion as illustrative the values 
|Am2^| = (2.0; 3.0) x 10"^ eV^ and sin^ 2^23 = 0.92; 1.0. 

^We have not written explicitly the two possible Majorana CP-violation phases [Ifil [T7j which do not 
enter into the expressions for the oscillation probabilities of interest jl6i 1I8] . 



Let us note that the atmospheric neutrino and K2K data do not allow one to determine 
the signs of ^m\^, and of cos 2^23 if sin^ 2^23 7^ l-O- This implies that in the case of 3- 
neutrino mixing one can have Amg^ > or Arrig^ < 0. The two possibilities correspond to 
two different types of neutrino mass spectrum: with normal hierarchy (NH), irti < 1712 < rris, 
and with inverted hierarchy (IH), m^ < rrii < 1712. The fact that the sign of cos 2^23 is not 
determined when sin^ 2^23 7^ 1-0 implies that when, e.g., sin^2^23 = 0.92, two values of 
sin^ 6^23 are possible, sin^6'23 = 0.64 or 0.36. 

The precise limit on the angle 9i^ from the CHOOZ and Palo Verde data is Am-g^ — 
dependent (see, e.g, ref. J22j). Using the 99.73% allowed range of Amg^ = (1.1 — 3.2) x 
10~^ eV^, from ref. |2S1, one gets from a combined 3-neutrino oscillation analysis of the solar 
neutrino, CHOOZ and KamLAND data ^: 

sin^^ia < 0.047 (0.074), 90% (99.73%) C.L. (5) 



The global analysis of the solar, atmospheric and reactor neutrino data performed in ref. 
gives sin^^is < 0.054 at 99.73% C.L. 

It is difficult to overestimate the importance of getting more precise information about 
the value of the mixing angle ^13, of determining the sign of Arrag^, or the type of the neutrino 
mass spectrum (with normal or inverted hierarchy), and of measuring the value of sin^ ^23 
with a higher precision, for the future progress in the studies of neutrino mixing. Although 
this has been widely recognized, let us repeat the arguments on which the statement is based. 

The mixing angle 6*13, or the absolute value of the element U^z of the PMNS matrix, 
\Uez\ = sin ^13, plays a very important role in the phenomenology of the 3-neutrino oscilla- 
tions. It drives the sub-dominant z/^ <->■ Ug {v^ <->■ z/g) oscillations of the atmospheric v^ {v^) 
and z/g (z/g) [SHI ED]- The value of Oi^ controls also the z/^ — »• z/g, t'/i — > z/g, z/g -^ v^ and z/g — > z/^ 
transitions in the long baseline neutrino oscillation experiments (MINOS, CNGS), and in 
the widely discussed very long baseline neutrino oscillation experiments at neutrino factories 
(see, e.g., refs. |^ EH ESI)- The magnitude of the T- violating and CP- violating terms in 
neutrino oscillations probabilities is directly proportional to sin ^13 (see, e.g., refs. jHHEHj)- 

If the neutrinos with definite mass are Majorana particles (see, e.g., ref. PI]), the pre- 
dicted value of the effective Majorana mass parameter in neutrinoless double /3— decay de- 
pends strongly in the case of normal hierarchical or partially hierarchical neutrino mass 
spectrum on the value of sin^6'i3 (see, e.g., ref. ISE])- 

The sign of ^fn^^ determines, for instance, which of the transitions (e.g., of atmospheric 
neutrinos) v^ — > z/g and z/g — *■ z/^, or u^ — *■ z/g and z/g -^ z/^, can be enhanced by the Earth 
matter effects |H7| EHl EH]- The predictions for the neutrino effective Majorana mass in 
neutrinoless double /3— decay depend critically on the type of the neutrino mass spectrum 
(normal or inverted hierarchical) [311 EH] • The knowledge of the value of 6'i3 and of the sign 
of Amg]^ is crucial for the searches for the correct theory of neutrino masses and mixing as 
well. 

Somewhat better limits on sin^ ^13 than the existing one can be obtained in the MINOS, 
OPERA and ICARUS experiments [111112] • Various options are being currently discussed 
(experiments with off-axis neutrino beams, more precise reactor antineutrino and long base- 
line experiments, etc., see, e.g., ref. [IS]) of how to improve by at least an order of magnitude, 
i.e., to values of ~ 0.005 or smaller, the sensitivity to sin^ 61^. The sign of ^Tn\^ can be de- 
termined in very long baseline neutrino oscillation experiments at neutrino factories (see. 



e.g., refs. |3H I32j). and, e.g, using combined data from long baseline oscillation experiments 
at the JHF facility and with off-axis neutrino beams |13] . If the neutrinos with definite mass 
are Majorana particles, it can be determined by measuring the effective neutrino Majorana 
mass in neutrinoless double /5— decay experiments |36| 1^ . Under certain rather special 
conditions it might be determined also in experiments with reactor z/g 45 . 

In the present article we study possibilities to obtain information on the value of sin^ ^13 
and on the sign of AttIj^ using the data on atmospheric neutrinos, which can be obtained in 
experiments with detectors able to measure the charge of the muon produced in the charged 
current (CC) reaction by atmospheric z/^ or z/^. It is a natural continuation of our similar 
study for water-Cerenkov detectors jlH]. In the experiments with muon charge identification 
it will be possible to distinguish between the z/^ and z/^ induced events. As is well known, 
the water-Cerenkov detectors do not have such a capability. Among the operating detectors, 
MINOS has muon charge identification capabilities for multi-GeV muons ^41|. The MINOS 
experiment is currently collecting atmospheric neutrino data. The detector has relatively 
small mass, but after 5 years of data-taking it is expected to collect about 440 atmospheric 
z/^ and about 260 atmospheric t'^ multi-GeV events (having the interaction vertex inside the 
detector). There are also plans to build a 30-50 kton magnetized tracking iron calorimeter 
detector in India within the India-based Neutrino Observatory (INO) project ^^. The 
INO detector will be based on MONOLITH design ^48^. The primary goal is to study the 
oscillations of the atmospheric u^j, and z/^. This detector is planned to have efficient muon 
charge identification, high muon energy resolution (~ 5%) and muon energy threshold of 
about 2 GeV. It will accumulate sufficiently high statistics of atmospheric z/^ and z/^ induced 
events in several years, which would permit to search for effects of the subdominant u^ — > z/g 
(z/g — > i^fj,) and z/^ -^ Ve {pe -^ i^fi) transitions. 

If Un and z/^ with energies E,y^p ^ 2 GeV take part in 2-neutrino u^ -^ Ur and u^ -^ i?r 
oscillations, one would have P{i^n -^ Vr) = P{i^fj. -^ ^r)- If a given detector collects a 
sample of atmospheric z/^ and z/^ induced CC events with /x^ having energies exceeding a 
few GeV, the difference between the /i~ and /i"^ event rates in magnetized iron detectors of 
interest (MINOS, INO), apart from detector effects, would essentially be due to the difference 
between the z/^ and z/^ charged current (CC) deep inelastic scattering cross sections. In 
the case of 3- neutrino oscillations of the atmospheric z/^, z/^, z/g and z/g, the Earth matter 
effects can generate a difference between the z/^ — *> z/g (z/g -^ z/^) and u^ -^ z/g (z/g -^ z/^) 
transitions. For Amlj^ -C |A?7i3^|, which is implied by the current solar and atmospheric 
neutrino data, and if sin^ ^13 7^ 0, the Earth matter effects can resonantly enhance either 
the z/^ — > z/g and z/g — > z/^, or the z/^ —>■ z/g and z/g — > z/^ transitions, depending on the sign 
of Am^i- The effects of the enhancement can be substantial for sin^ ^13 ^ 0.01. Under 
the condition Am^i ^ jAmgJ, the expressions for the probabilities of z/^ —^ z/g (z/g -^ z/^) 
and z/^ — > z/g (z/g —>■ v^) transitions contain also sin^ 6*23 as a factor which determines their 
maximal values. Since the fluxes of multi-GeV z/^ and z/g (z/^ and z/g) differ considerably 
(by a factor ~ (2.6 — 4.5) for E^^y ~ (2 — 10) GeV), the Earth matter effects can create a 
substantial difference between the /i~ and /i"*" rates of events, produced by the atmospheric 
z/^ and z/^ in MINOS, INO, or any other detector of the same type. This difference can be 
relatively large and observable in the samples of the multi-GeV /i^ events (i?^ ~ (2 — 10) 
GeV), in which the muons are produced by atmospheric neutrinos with relatively large path 
length in the Earth, i.e., by neutrinos crossing deeply the Earth mantle [S21ISS] or the mantle 



and the core [23 EOl ED E2] • These are the neutrinos for which jlSl cos^„ ^ (0.3 — 0.4), 9n 
being the Nadir angle characterizing the (atmospheric) neutrino trajectory in the Earth. 

For Ar77,|^ > 0, the u^ —>■ v^ {p^ — > v^ and z/g — > z/^ {pe -^ ^^ transitions in the Earth 
lead to a reduction of the rate of the multi-GeV \i~ events observable in MINOS, INO, etc., 
with respect to the case of absence of these transitions (see, e.g., refs. [^ W^ I3U1 1^ 13^). 
If Arrig^ < 0, the [i^ event rate will be reduced. Correspondingly, as observable which is 
sensitive to the Earth matter effects, and thus to the value of sin^ d\j, and the sign of Amg^^, as 
well as to sin^ 6*23, we can consider the Nadir-angle distribution of the ratio A^(/i^)/A^(/i^) of 
the multi-GeV \i~ and /x"*" event rates, or, equivalently, of the [i~ — /i"*" event rate asymmetry 
^M"M+ ~ [^(a^~) ~ ^(a*^)]/[^(/^ ) + ^(a^"^)]- The systematic uncertainty, in particular, in 
the Nadir angle dependence of N{fi~)/N{^'^), and correspondingly in the asymmetry A^-^+, 
can be smaller than those on the measured Nadir angle distributions of the rates of /i~ and 
/i"*" events, N{fj,~) and N^fi'^). 

We have obtained predictions for the Nadir-angle distribution of ^4^-^+ in the case of 
3- neutrino oscillations of the atmospheric z/^, z/^, z/g and z/g, both for neutrino mass spectra 
with normal {Aml^ > 0) and inverted {Aml^ < 0) hierarchy, (A^'l +)nh5 (^u-u+)ih! and for 
sin^023 = 0.64; 0.50; 0.36. These are compared with the predicted Nadir-angle distributions 
of the same asymmetry in the case of 2-neutrino z/^ —>■ Vr and z/^ — > Vr oscillations of 
the atmospheric z/^ and z/^ (i.e., for sin^ 6*13 = 0), Aj^- +. Predictions for the three types 
of asymmetries indicated above of the suitably integrated Nadir angle distributions of the 
/i~ and /i"*" multi-GeV event rates are also given. Our results show, in particular, that for 
sin^ 6*23 ^ 0.50 and sin^ 26'i3 ^ 0.06 the effects of the Earth matter enhanced subdominant 
transitions of the atmospheric neutrinos, u^ — > z/g and and z/g — i> z/^, or u^ -^ z/g and z/g — ;> z/^, 
can be sufficiently large to be observable with INO and possibly with MINOS detectors. 
Conversely, if the indicated effects are observed in the MINOS and/or INO experiments, 
that would imply that sin^ 26'i3 ^ 0.05, sin^ ^23 ~ 0.50 and at the same time would permit 
to determine the sign of Am|^ and thus to answer the fundamental question about the type 
of hierarchy - normal or inverted, the neutrino mass spectrum has. 

Let us note that the Earth matter effects in atmospheric neutrino oscillations have been 
widely studied (for recent detailed analyzes see, e.g., refs. [221 HE], which contain also a 
rather complete list of references to earlier work on the subject). A rather detailed analysis 
for the MONOLITH detector has been performed in ref. ^^. A large number of studies 
have been done for the Super-Kamiokande detector, or more generally, for water-Cerenkov 
detectors. In ref. jlH] , in particular, the magnitude of the Earth matter effects in the Nadir 
angle distribution of the ratio of the multi-GeV /i— like and e— like events, measured in 
water-Cerenkov detectors, N^/Nf., has been investigated. This Nadir angle distribution is 
the observable most sensitive to the matter effects of interest. It was concluded that for 
sin^ 6*23 ^ 0.50, sin^ ^13 ^ 0.01 and Amg^ > 0, the effects of the Earth matter enhanced 
T^fi -^ ^e and z/g -^ Vy, transitions of the atmospheric v^ and z/g, might be observable with 
the Super-Kamiokande detector. However, determining the sign of Am|^ would be quite 
challenging in this experiment (or its bigger version - Hyper-Kamiokande 54 ). In general, 
the matter effects in the Nadir angle distribution of the ratio of the multi-GeV /i— like and 
e— like events, N^/Nf,., which can be measured in the Super-Kamiokande or other water- 
Cerenkov detectors, are smaller than the matter effects in the Nadir angle distribution of the 
ratio of the multi-GeV /x~ and /x"*" events, N{fj,~)/N{fj,~^), which can be measured in MINOS, 



INO, or any other atmospheric neutrino experiment with sufficiently good muon charge 
identification. The reason is that in the case of water-Cerenkov detectors, approximately 
2/3 of the rate of the multi-GeV /x — like events is due to p^, and ~ 1/3 is due to i/^; 
similar partition is valid for the multi-GeV e— like events. Depending on the sign of Am|]^, 
the matter effects enhance either the neutrino transitions, z/^ —>■ v^ and and z/g — > z/^, or 
the antineutrino transitions, v^i_ —>■ z/g and Ue —>■ i^fi, but not both types of transitions. 
Correspondingly, because the z/^,e— and z/^,e— induced events are indistinguishable in water- 
Cerenkov detectors, only ~ 2/3 or ~ 1/3 of the events in the multi-GeV yU— like and e— like 
samples collected in these detectors are due to neutrinos whose transitions can be enhanced 
by matter effects. This effectively reduces the magnitude of the matter effects in the samples 
of multi-GeV /i— like and e— like events. Obviously, such a "dilution" of the magnitude of 
the matter effects does not take place in the samples of the multi-GeV fi~ and /x"*" events, 
which can be collected in MINOS and INO experiments, i.e., in the experiments with muon 
charge identification. 

2 Subdominant 3-z/ Oscillations of Multi-GeV Atmo- 
spheric Neutrinos in the Earth 

In the present Section we review the physics of the subdominant 3-neutrino oscillations 
of the multi-GeV atmospheric neutrinos in the Earth (see, e.g., ref. |46j ) . 

The subdominant z/^ -^ z/g (z/^ -^ z/g) and z/g -^ z/^(^) (z/g -^ ^^i{t)) oscillations of the multi- 
GeV atmospheric neutrinos of interest should exist and their effects could be observable if 
three-flavor-neutrino mixing takes place in vacuum, i.e., if sin^ 2^13 7^ 0, and if sin^ 2^13 is 
sufficiently large |2nil^E21- These transitions are driven by /\m\-^. The probabilities of 
the transitions contain sin^ 62^ as factor which determines their maximal value (see further) . 
For Amg]^ > 0, the v^ — > z/g (z/^ -^ z/g) and z/g -^ ^f^ir) (j^e -^ ^K^)) transitions of the 
multi-GeV atmospheric neutrinos, are amplifled (suppressed) by the Earth matter effects; 
if Am|^ < 0, the transitions of neutrinos are suppressed and those of antineutrinos are 
enhanced. Therefore for a given sign of Am|^, the Earth matter affects differently the 
transitions of neutrinos and antineutrinos. Thus, the study of the subdominant atmospheric 
neutrino oscillations can provide information, in particular, about the type of neutrino mass 
hierarchy (or sign of Ar/ig^), the magnitude of sin^ ^13 and the octant where 623 lies. 

Under the condition lArrig^l ^ ^fn^i, which the neutrino mass squared differences de- 
termined from the existing atmospheric and solar neutrino and KamLAND data satisfy, the 
relevant 3-neutrino z/^ -^ z/g (z/^ -^ z/g) and z/g -^ u^^r) ipe — ^ ^m(t)) transition probabilities 
reduce effectively to a 2-neutrino transition probability [55J with Am|^ and ^13 playing the 
role of the relevant 2-neutrino oscillation parameters. 

The 3-neutrino oscillation probabilities of interest for atmospheric z/g^^ having energy 
E and crossing the Earth along a trajectory characterized by a Nadir angle On-, have the 
following form [S3j: 

Pzui^T^e ^ Z/g) = 1 - P2., (6) 

P^M ^ y^.) = P3u{l^^ -^ ^e) = Sl^ P2u, (7) 

P3u{l^e ^ l^r) = 4, P2,, (8) 



PU^^. ^ z/^) = 1 - 4 P2. - 24s^3 [1 - Re {e'^^A^M ^ i^r))] , (9) 

Here P21, = P2i/(^^^3i' ^i35 -^' ^n) is the probability of 2-neutrino z/g — > z/^ oscillations in the 
Earth, where z^^ = 8231^^ + C231^t 55 , and k and A2u{j^t ~^ ^r) = ^2^ are known phase and 
2-neutrino transition probability amplitude [23 EDI ESI- 

The fluxes of atmospheric z/g,^ of energy E, which reach the detector after crossing the 
Earth along a given trajectory specified by the value of 6*^, $y^^(£^, 6'„), are given by the 
following expressions in the case of the 3-neutrino oscillations under discussion pUj : 



^,XE.er.) = <^l {l + [sl^r-l]P2.), (11) 

$,^(E, On) = $° (1 + 4 [(^23 r)-' - 1] P2U - 2cl,sl, \l - Re {e-^^A2M ^ ^r)) 



(12) 
where $° = $° {E, On) is the Veip) flux in the absence of neutrino oscillations and 

The interpretation of the SK atmospheric neutrino data in terms of z/^ -^ v^ oscillations 
requires the parameter S23 to lie approximately in the interval (0.30 - 0.70), with 0.5 being 
the statistically preferred value. For the predicted ratio r{E, On) of the atmospheric z/^ 
and Ve fluxes for i) the Earth core crossing and ii) only mantle crossing neutrinos, having 
trajectories for which 0.3 < cos^„ < 1.0, one has fl^EZlEH]: r{E,0^) ^ (2.0 - 2.5) for 
the neutrinos giving contribution to the sub-GeV samples of Super-Kamiokande events, and 
r{E,On) = (2.6 — 4.5) for those giving the main contribution to the multi-GeV samples. If 
5^3 = 0.5 and r{E,0^) ^ 2.0, we have (5^3 r{E,0,) - 1) ^ 0, ((5^3 r{E,0^))-^ - 1) ^ 0, 
and the possible effects of the z/^ — > u^ and z/g —>■ v^ir) transitions on the z/g and z/^ fluxes, 
and correspondingly in the sub-GeV e— like and /x— like samples of events, would be strongly 
suppressed independently of the values of the corresponding transition probabilities. For 
the multi-GeV neutrinos one flnds s^r^ll - (s^g r{E,0,))-^] = 0.06 - 0.14 (0.16 - 0.27) and 
(s^g r{E,0,) - 1) = 0.3 - 1.3 (0.66 - 1.9) for 5^3 = 0.5 (0.64). Obviously, the effects of 
interest are much larger for the multi-GeV neutrinos than for the sub-GeV neutrinos. They 
are also predicted to be larger for the flux of (and event rate due to) multi-GeV atmospheric 
z/g than for the flux of (and event rate due to) multi-GeV atmospheric z/^. 

The same conclusions are valid for the effects of oscillations on the fluxes of, and event 
rates due to, atmospheric antineutrinos z/g and z/^. The formulae for anti-neutrino fluxes and 
oscillation probabilities are analogous to those for neutrinos: they can be obtained formally 
from eqs. (jHl) - (fTT^ by replacing the neutrino related quantities - probabilities, k, A2u{i^t -^ 
Vt) and fluxes, with the corresponding quantities for antineutrinos: P2v{^fn\i, O13] E, 0„) —>■ 

P2^{Aml^,0i3;E,0n), K -^ R, A2u{i^t -^ ^t) -^ A2„{i^r -^ i^r) = Au, P3u{^l "^ J^l') ^ 

P3M ^ i^i'), ^^^liE,On) -^ <^fl{E,On) and r(E,^„,) -^ f (E, 0„,) (see refs. 130111^1). 

Equations © - ©7 (fTT|) - (fT^ and the similar equations for antineutrinos imply that 
in the case under study the effects of the z/^ — *> z/g, z/^ -^ z/g, and z/g — >■ z/;i(r), t'e -^ '^m(t); 
oscillations i) increase with the increase of S23 and are maximal for the largest allowed value 
of 5^3, ii) should be substantially larger in the multi-GeV samples of events than in the 



sub-GeV samples, and iii) in the case of the multi-GeV samples, for Ar7i|^ > they lead to 
a decrease of the fi~ event rate, while if AttIj^ < 0, the /x"'" event rate will decrease. The last 
point follows from the fact that the magnitude of the effects we are interested in depends 
also on the 2- neutrino oscillation probabilities, P2u and P2u, and that P2u or P2u (but not 
both probabilities) can be strongly enhanced by the Earth matter effects. In the case of 
oscillations in vacuum we have P2u = P2u ~ sin^ 2^13. Given the existing limits on sin^ 26'i3, 
the probabilities P2U and P2U cannot be large if the oscillations take place in vacuum. 

If sin^ ^13 7^ 0, the Earth matter effects can resonantly enhance either the u^ — > Ue and 
^e -^ ^fi, or the Ufj, -^ z/g and z/g — > t'^ transitions, depending on the sign of Am^^. The 
enhancement mechanisms are discussed briefly in the next subsection. 

2.1 Enhancing Mechanisms 

As is well-known, the Earth density distribution in the existing Earth models is assumed 
to be spherically symmetric ^ and there are two major density structures - the core and the 
mantle, and a certain number of substructures (shells or layers). The core radius and the 
depth of the mantle are known with a rather good precision and these data are incorporated in 
the Earth models. According to the Stacey 1977 and the more recent PREM models [5^l6Uj. 
which are widely used in the calculations of the probabilities of neutrino oscillations in the 
Earth, the core has a radius Re = 3485.7 km, the Earth mantle depth is approximately 
Rman = 2885.3 km, and the Earth radius is R^ = 6371 km. The mean values of the 
matter densities and the electron fraction numbers in the mantle and in the core read, 
respectively: p^an = 4.5 g/cm^, p^ ^ 11.5 g/cm^, and [HI] i;*""" = 0.49, y/ = 0.467. 
The corresponding mean electron number densities in the mantle and in the core read: 
N^""^ = pmanYJ^''''/mN = 2.2 Nactti-^ , N^ = PcYe/mN = 5.4 NAcm-^, rriN and Na being 
the nucleon mass and Avogadro number. 

Numerical calculations show [23 El] that, e.g., the z/g — »• z/^ oscillation probability of in- 
terest, calculated within the two-layer model of the Earth with pman (or iV^"") and pc (or N^) 
for a given neutrino trajectory determined using the PREM (or the Stacey) model, repro- 
duces with a remarkably high precision the corresponding probability, calculated by solving 
numerically the relevant system of evolution equations with the much more sophisticated 
Earth density profile of the PREM (or Stacey) model. 

In the two-layer model, the oscillations of atmospheric neutrinos crossing only the Earth 
mantle (but not the Earth core), correspond to oscillations in matter with constant density. 
The relevant expressions for P2u, k and A2y(z/,- — > z/,-) are given by (see, e.g., ref. |l?)j): 

P2u{Aml, 013; E, On) = sin^ (^^\ sin^ 2^, (14) 

-^\[^L + V-2G,N—L - ^], (15) 

A2,{Ur ^ yr) = I + fe-'^^-l^COS^^;, (16) 



■^Let us note that because of the approximate spherical symmetry of the Earth, a given neutrino trajectory 
through the Earth is completely specified by its Nadir angle. 
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Here 

\ 2 



AM^ = Am|i 



\ 



1-^1 cos= 26,3 + sill' 26,3 , (17) 

rman I 



is the mass difference between tlie two mass-eigenstate neutrinos in tfie mantle, Q'^ is tfie 
mixing angle in the mantle, 



■2ryni _ sin' 2^13 

sm W^ 



;i_£^)2cos2 2^13 + sin^ 2^13 

L is the distance the neutrino travels in the mantle, pman and p^^a^ ^^^ the mean density 
along the neutrino trajectory and the resonance density in the mantle, 

_ Am^i cos 2^13 

For a neutrino trajectory which is specified by a given Nadir angle Qn we have: 

L = 2i?® cos On (20) 

where R^ = 6371 km is the Earth radius (in the PREM [00] and Stacey j^HI models) ^. 

Consider for definiteness the case of Am|^ > 0. It follows from eqs. pij) and (fT^ that 
the oscillation effects of interest will be maximal if P2u — 1- The latter is possible provided 
i) the well-known resonance condition [^EHI, leading to sin' 2^^ = 1, is fulfilled, and ii) 
cos (^If^) = -1. Given the values of pman and F™'*", or the value of iV™"", the first 
condition determines the neutrino energy at which P21, can be enhanced: 

E„. ^ 6.6 (,^^] (^^^] cos 29,3 GeV . (21) 



man 



10-3 eV2; \^ N 

If the first condition is satisfied, the second determines the length of the path of the neutrinos 
in the mantle for which one can have P2u — 1: 

= 1.27r tan 2^13 ^ -]{^— — = tt, 22 

V 2E J^,^^ \NAcm-^J Uo^km; ^ ^ 

Taking Amj-^ ^ (2.0 - 3.0) x lO'^ eV, N^^'^ ^ 2 NAcm'^ and cos 2^13 ^ 1 one finds from 
eq. (J2H): Eres — (6.6 — 10.0) GeV. The width of the resonance in E, 26E, is determined, as 
is well-known, by tan 2^13: 6E/Eres ~ tan26'i3. For sin' ^13 ~ (0.01 — 0.05), the resonance 
is relatively wide in the neutrino energy: 5E / E^es — (0.27 — 0.40). Equation (j221) implies 
that for sin' ^13 = 0.05 (0.025) and N^'''' = 2.2 NAcm^^ one can have Pgz. = 1 only if 
L ^ 8000 (10000) km. 

It follows from the above simple analysis [33J that the Earth matter effects can amplify 
P2iy significantly when the neutrinos cross only the mantle i) for E ^ (6 — 11) GeV, i.e., in 
the multi-GeV range of neutrino energies, and ii) only for sufficiently long neutrino paths in 

^Neutrinos cross only the Earth mantle on the way to the detector if 9n ^ 33.17°. 
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the mantle, i.e., for cos^„ ^ 0.3. The magnitude of the matter effects of interest increases 
with increasing of sin^ ^13. 

In the case of atmospheric neutrinos crossing the Earth core, new resonant effects become 
apparent. For sin^^is < 0.05 and Am|^ > 0, we can have P21/ — 1 only due to the effect of 
maximal constructive interference between the amplitudes of the the z/g — ^ u'^ transitions in 
the Earth mantle and in the Earth core 123 UHl HHI ■ The effect differs from the MSW one j^H] 
and the enhancement happens in the case of interest at a value of the energy between the 
resonance energies corresponding to the density in the mantle and that of the core. The 
mantle-core enhancement effect is caused by the existence (for a given neutrino trajectory 
through the Earth core) oi points of resonance-like total neutrino conversion, P21, = 1, in the 
corresponding space of neutrino oscillation parameters I^EHI- The points where P21/ = 1 
are determined by the conditions [13 EO] : 



tan0-±J 77'^-^ „ tan0- = ± , '''^^- (23) 

Vcos(2^^-4^J' '^ ./-cos2^;;,cos(2^;;-4^J ^ ^ 



where the signs are correlated and cos26'^ cos (26'^— 4^^) < 0. In eq. pHjl 2(f)' and 2(f)" are the 
oscillation phases (phase differences) accumulated by the (two) neutrino states after crossing 
respectively the first mantle layer and the core, and 6'^ is the neutrino mixing angle in the 
core. A rather complete set of values of Am'^^/E and sin^ 26'i3 for which both conditions 
in eq. ()23p hold and P2U = 1 for the Earth core-crossing atmospheric z/^ and z/g having 
trajectories with Nadir angle 6n = 0; 13°; 23*^; 30*^ was found in ref. jHOj- The location of 
these points determines the regions where P2U is large, P2U ~ 0.5. These regions vary slowly 
with the Nadir angle, they are remarkably wide in the Nadir angle and are rather wide in 
the neutrino energy jSU], so that the transitions of interest produce noticeable effects: we 
have 6E/E = 0.3 for the values of sin^ 6*13 of interest jSOl EHI- 

For sin^ ^13 < 0.05, there are two sets of values of Ar7i|]^ and sin^ 613 for which eq. pH|l 
is fulfilled and P2U = 1- These two solutions of eq. (j^ occur for, e.g., 9n = 0; 13°; 23°, 
at 1) sin2 2^i3 = 0.034; 0.039; 0.051, AmlJE = 7.2; 7.0; 6.5 x 10"^ eVVMeV, and at 
2) sin2 20i3 = 0.15; 0.17; 0.22, AmlJE = 4.8; 4.5; 3.8 x 10"^ eVVMeV (see Table 2 
in ref. jSD])- The first solution corresponds to jSH] cos 20' = —1, cos 20" = — 1 and ^ 
sin^(26'^ - 46'^) = 1. For Am^^ = 2.0 (3.0) x 10"^ eV^, the total neutrino conversion 
occurs in the case of the first solution at E ^ (2.8 - 3.1) GeV {E = (4.2 - 4.7) GeV). The 
values of sin^ 2^13 at which the second solution takes place are marginally allowed. If, e.g., 
ATTig^ = 2.5 X 10~^ eV^, one has P2U = 1 for this solution for a given 6n in the interval 
< ^„ < 23° at E lying in the interval E = (5.3 - 6.7) GeV. 

The effects of the mantle-core enhancement of P21/ (or P2u) increase rapidly with sin^ 26'i3 
as long as sin^ 26'i3 ^ 0.06, and should exhibit a rather weak dependence on sin^ 2^13 for 
0.06 ^ sin^ 26'i3 < 0.19. If 3-neutrino oscillations of atmospheric neutrinos take place, 
the magnitude of the matter effects in the multi-GeV /i— like and e— like event samples, 
produced by neutrinos crossing the Earth core, should be larger than in the event samples 
due to neutrinos crossing only the Earth mantle (but not the core). This is a consequence of 
the fact that in the energy range of interest the atmospheric neutrino fluxes decrease rather 



^The term "neutrino oscillation length resonance" (NOLR) was used in ref. [IHI to denote the mantle-core 
enhancement effect in this case. 
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rapidly with energy - approximately as E~'^-'^, while the neutrino interaction cross section 
rises only linearly with E, and that the maximum of P2U (or P2iy) due to the resonance- 
like mantle-core interference effect takes place at approximately two times smaller energies 
than that due to the MSW effect for neutrinos crossing only the Earth mantle (e.g., at 
E ^ (3.5 - 3.9) GeV and E = 8.3 GeV, respectively, for Am|i = 2.5 x IQ-^ eV^). 

The same results, eqs. ()21|) and ()22|) . and conclusions are valid for the antineutrino oscil- 
lation probability P2U in the case of Amg^^ < 0. As a consequence, a preferable detector for 
distinguishing the type of mass hierarchy would be the one with muon charge discrimination, 
such that neutrino interactions can be distinguished from those due to antineutrinos. 

3 Results 

It follows from the preceding analysis that in the case of detectors with muon charge 
identification, as observable which is most sensitive to the Earth matter effects, and thus 
to the value of sin^ 6*13 and the sign of Amli, as well as to sin^6'23, we can consider the 
Nadir-angle (6'„) distribution of the ratio iV(/x~)/A^(/i+) of the multi-GeV fi" and /x"*" event 
rates, or equivalently the Nadir-angle distribution of the /i~ — /i"*" event rate asymmetry 



Njfi-) - iV(^^ 
N{^i-) + N{fi^ 



A,-,^ = zr.i 7 ::rj ■ m 



We have obtained predictions for the cos6'„ distribution of the ratio N{fi~)/N{fi~^) and the 
asymmetry ^^-/x+ in the case of 3-neutrino oscillations of the atmospheric z/^, z/^, u^ and z/g, 
both for neutrino mass spectra with normal (Am|^ > 0) and inverted (Am|]^ < 0) hierarchy, 
and for sin^ 6^23 = 0.64; 0.50; 0.36, and sin^ 26^13 = 0.05; 0.10. These are compared with 
the predicted Nadir-angle distributions of the same ratio and asymmetry in the case of 2- 
neutrino (sin^ ^13 = 0) vacuum u^ -^ u^ and z/^ -^ Vr oscillations of the atmospheric z/^ and 

In the calculations we have used the predictions for the Nadir angle and energy distribu- 
tions of the atmospheric neutrino fluxes given in ref . |2H1 • The interactions of the atmospheric 
neutrinos are described by taking into account only the z/^ and v^ deep inelastic scattering 
(DIS) cross sections. The latter are calculated using the GRV94 parton distributions given 
in ref. [021 • We present here results for the asymmetry ^/x-^+ ^ ■ They are shown graph- 
ically in Figs. H - ini The figures correspond to three different intervals of integration over 
the energies of the atmospheric v^ and z/^, and of the /i~ and [i^ they produce in the de- 
tector, E = [2, 10], [2, 20], [5, 20] GeV, and ^ thus to three different possible event samples. 
Figures HI - m show the the asymmetries A^~^+ and ^^-^+ as functions of cos6'„ for two 
"reference" values of |Am|i|, Am^^ = ±2 x 10~^ eV^ and Aml^ = ±3 x 10~^ eV^, while in 

^It is interesting to note that the ratio Af(/i^)/A^(/i+) exhibits essentially the same dependence on 
cos^n as the asymmetry A^-^+. This is a consequence of the special form of the dependence of A^-^+ 
on N{iJ.^)/N{fi^) and of the fact that typically one finds N{fi^)/N{jj,^) ^ (1.5 — 2.4) for the ranges of 
the values of the parameters of interest. Correspondingly, the following approximate relation holds (within 
~ 20% and typically with much higher precision) for the range of values of the parameters of interest: 

^The iron-magnetized calorimeter detectors allow to reconstruct with a certain precision the initial neu- 
trino energy as well, see refs. 41 , 47j, :48.; . 
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Figs. El-Elwe present results for the asymmetries in the rates of the multi-GeV fi~ and /i"^ 
events, integrated over cos6'„ in the intervals [0.30,0.84] {mantle bin) and [0.84,10] {core bin), 
A^-^+ and ^^-^+- The dependence of the latter on sin^ 2^13 for Amg^ = ±2 x 10~^ eV^, 
and on Aml^ for sin^ 26'i3 = 0.10, is shown for three values of sin^^23 = 0.36; 0.50; 0.64. 

As Figs. ^ - ini indicate, the Earth matter effects can produce a noticeable deviations 
of A^-^+ from the 2-neutrino vacuum oscillation asymmetry ^^-„+ at cos0„ ^ 0.3. As a 
quantitative measure of the magnitude of the matter effects one can use the deviation of the 
asymmetry Afj_-^+ in the case of 3-neutrino oscillations, sin^ 26'i3 7^ 0, sin^ 2^13 ^ 0.04, from 
the asymmetry, ^4^1+, predicted in the case of 2-neutrino oscillations, i.e., for sin^ 26'i3 = 0, 
or the relative difference between the two asymmetries, 

A = ^---:,, ^---^ . (25) 

The magnitude of the matter effects, or the relative difference A, depends critically on the 
value of sin^ ^23 : | A | increases rapidly with the increasing of sin^ ^23 • This is clearly seen 
in Figs. H-|ni The matter effects in A^-^+ are hardly observable for sin^ ^23 ~ 0.30. For 
cos6'„ < 0.84, i.e., in the mantle bin, the asymmetry difference |A| increases practically 
linearly with sin^2^i3. In the case of 0.84 < cos^„ < 1.0, i.e., in the core bin, |A| increases 
rapidly with sin^2^i3 until the latter reaches the value of sin^ 2^13 = 0.06. For values of 
sin^ 2^13 = (0.06 — 0.15), A is essentially independent of sin^ 2^13 and is given by its value at 
sin^ 26'i3 — 0.06 (Figs. El- EI)- The magnitude of the asymmetry difference A depends weakly 
on Am^i taking values in the interval ~ (2 — 3) x 10"'^ eV^, as long as the energy integration 
interval is sufficiently wide to include the energy regions where the Earth matter effects 
enhance strongly the subdominant transition probabilities. If this is not the case, a noticeable 
dependence on A?ti|^ can be present. This is illustrated e.g, in Fig. which corresponds to 
E = [2, 10] GeV. The asymmetry difference in the mantle bin diminishes monotonically as 
|Am|;^| increases starting from the value of ~ 1.3 x 10~^ eV^ and becomes rather small at 
lAmg;^! ^ 3 X 10^^ eV^. This behavior can be easily understood: for |A?7i3j^| = 2 x 10^^ eV^, 
the region of enhancement of the subdominant neutrino oscillations lies in the region of 
energy integration, while for |Am|]^| > 3 x 10~^ eV^ the enhancement region is practically 
outside the region of integration over the neutrino energy. 

For the ranges considered of the three oscillation parameters, sin^ 6'23, sin^ 2^13 and 
|Am|]^|, the magnitude of the asymmetry difference |A| depends weakly on the maximal 
neutrino (and muon) energy, Emax, for the chosen event sample as long as E^ax ~ 10 GeV. 
By increasing the minimal energy of the neutrinos contributing to the event sample, Emm, 
from 2 GeV to, e.g., 5 GeV, one could diminish the asymmetry in the core bin substantially 
(Fig. 13]). In that case, a large fraction of the region of enhancement is not included within 
the interval of integration. 

For sin2^23 ^ 0.50, sin^ 2^13 > 0.06 and |Am|i| = (2 - 3) x lO^^ gv2, the Earth matter 
effects produce an integrated asymmetry difference |A| which is bigger than approximately 
~ 15%, can reach the values of (30 — 35)% (Figs. El -El), and thus can be sufficiently large to 
be observable. As Figs. El- El clearly show, the sign of the relative difference of the integrated 
asymmetries is anticorrelated with the sign of Arn^^: for Am^i > we have A^-^+ < A^^- +, 
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while for A?7i|^ < 0, the inequahty Afj^-^+ > ^^-„+ holds ®. Therefore the measurement 
of A^-fj_+ can provide a direct information on the sign of Ar/ig^, i.e., on the neutrino mass 
hierarchy. 

It follows from Figs.^|3]that the deviations of the asymmetry ^^-/x+ from the 2-neutrino 
oscillation one, A^l +, are maximal typically i) in the "core bin", cos6'„ = [0.84 — 1.0], and 
ii) at cos6'„ ~ 0.50 in the "mantle bin", cos6'„ < 0.84. For Am^^ = 2 x lO^^ eV^ sin^ 6^23 = 
0.50 (0.64) and E = [2,10] GeV (Fig. d), we have at cos^„ ~ 0.50 and sin^ 2^13 = 0.05: 
Af_^+ = 0.29, while the matter effects lead to A^-^+ = 0.25 (0.24). This corresponds to a 
negative relative difference between y4(/i~/i"'") and A^''(/i~/i"'"), A < 0, and |A| ~ 14% (17%). 
For sin^ 26'i3 = 0.10, one finds A ~ —28% (34%). The relative difference A in the core 
bin is also negative and |A| has similar or larger values. For E in the interval E = [5,20] 
GeV (Fig. HI), we get at cos^„ ~ 0.50 for Amj^ = 3 x 10"^ eV^ sin^ 2^13 = 0.10 and 
sin^ ^23 = 0.64: A?^-n+ — 0.25, ^^-/x+ — 0.075, and a relative difference between the two 
asymmetries A = —70%. 

Our results show that the Earth matter effects in the Nadir-angle distribution of the ratio 
N{fi~)/N{fi~^) of the rates of multi-GeV fi~ and /x"*" events, or equivalently in the Nadir-angle 
distribution of the /i~ — /i"*" event rate asymmetry A^-^+, eq. (J2l|l . can be sufficiently large 
to be observable in the current and planned experiments with iron magnetized calorimeter 
detectors which have muon charge identification capabihties (MINOS, INO, etc.). 



4 Conclusions 

We have studied the possibilities to obtain information on the values of sin^ ^13 and 
sin^^23, and on the sign of Am|]^ using the data on atmospheric neutrinos, which can be 
obtained in experiments with detectors able to measure the charge of the muon produced 
in the charged current (CC) reaction by atmospheric u^ or z/^ (MINOS, INO, etc). The 
indicated oscillation parameters control the magnitude of the Earth matter effects in the 
subdominant oscillations, z/^ -^ Ve {^e -^ ^fi) and u^ — > z/g (z/g -^ i^p.), of the multi-GeV 
(i? ~ (2 — 10) GeV) atmospheric neutrinos. As observable which is most sensitive to the 
Earth matter effects, and thus to the value of sin^ ^13 and and the sign of Am|]^, as well as to 
sin^6'23, we have considered the Nadir-angle (^„) distribution of the ratio iV(/i~)/A^(/i+) of 
the multi-GeV fi~ and /x"*" event rates, and the corresponding fi~ — [i'^ event rate asymmetry 
^/i-/i+) 6q. ((211) • The systematic uncertainty, in particular, in the Nadir angle dependence 
of N{ii~)/N{ii^) and of the asymmetry A^-^+, can be smaller than those on the measured 
Nadir angle distributions of the rates of /i~ and /i"*" events, N{fi~) and N{fi~^). We have 
obtained predictions for the cos6'„ distribution of the asymmetry A^~^+ (and of the ratio 
N{fi~)/N{fj,~^)) in the case of 3-neutrino oscillations of the atmospheric z/^, u^, v^ and z/g, 
both for neutrino mass spectra with normal (AttIj^ > 0) and inverted (Arrigj^ < 0) hierarchy, 
and for sin^ ^23 = 0.64; 0.50; 0.36, and sin^ 2^13 = 0.05; 0.10. These are compared with the 
predicted Nadir-angle distribution of the same ratio and asymmetry in the case of 2-neutrino 
(sin^ ^13 = 0) vacuum z/^ -^ v^ and z/^ -^ Vr oscillations of the atmospheric v^ and z/^, A'^'i. +. 



^'Note that A _ + > 0. This is a consequence of the fact that the j/^ DIS cross section is approximately 
by a factor 2 bigger than the P^ DIS cross section and that the fluxes of atmospheric i^^ and P^, {i^e and i/g) 
do not differ considerably. 
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Our results are summarized in Figs. ^ - IHl Figures ^ - |3] show the dependence of the 
asymmetries Afj^-^+ and ^^-„+ on cos6'„ for two "reference" values of jAmg^l, Amg^ = 
±2 X 10~^ eV^ and Amgj^ = ±3 x 10~^ eV^, and for three possible ranges of energies of the 
atmospheric neutrinos, contributing to the event rates of interest, E = [2, 10], [2, 20], [5, 20] 
GeV. In Figs. El-Elwe present results for the asymmetries in the rates of the multi-GeV /i~ 
and n'^ events, integrated over cos 6'n in the intervals [0.30,0.84] {mantle bin) and [0.84,10] 
{core bin). We find that for sin^ ^23 ^ 0.50, sin^ 2^13 > 0.06 and jAm^J = (2-3) x lO'^ eV^, 
the Earth matter effects produce a relative difference between the integrated asymmetries 
^fi-fi+ ^-iid ■^'u-u+ which is bigger in absolute value than approximately ~ 15%, can reach the 
values of (30 — 35)% (Figs. El -El), and thus can be sufficiently large to be observable. As our 
results show (Figs. El -El), the sign of the indicated asymmetry difference, (A^-^+ — Aj^- +), 
is directly related to the sign of Amg^: for Aml^ > we have (y4^-^+ — A^l +) < 0, while 
if Am^i < then (^4^-^+ — ^4^1+) > Therefore the measurement of the Nadir angle 
dependence of ^4^-^+, or of the value of ^/i-/i+ in the mantle and/or in the core bins, can 
provide a direct information on the sign of A?7i|^, i.e., on the neutrino mass hierarchy. 

To summarize, the studies of the oscillations of the multi-GeV atmospheric u^^ and z/^ 
in experiments with detectors having good muon charge identification capabilities (MINOS, 
INO, etc.), can provide fundamental information on the values of sin^^is and sin^6'23, and 
on the sign of Aml^, i.e., on the neutrino mass hierarchy. 
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Figure 1: The Nadir angle distribution of the charge asymmetry, A^-^+, eq. (j23j), of the 
multi-GeV [i~ and /i"*" event rates, integrated over the neutrino (and muon) energy in the 
interval E = (2.0 — 10.0) GeV, in the cases i) of 2-neutrino u^ -^ Vr and z/^ — *> Vr os- 
cillations in vacuum of the atmospheric u^ and v^j_ and no z/e and v^ oscillations, ^^-o+ 
(solid lines), ii) 3- neutrino oscillations of i/^, i/^ v^ and z/g in the Earth and neutrino 
mass spectrum with normal hierarchy (Aj^'i +)nh (dashed lines), or with inverted hier- 
archy, (^i^-„+)iH (dotted lines). The results shown are for |Am3;^| = 2 x 10~^ eV^, 
sin^ 6*23 = 0.36 (upper panels); 0.50 (middle panels); 0.64 (lower panels), and sin^ 26'i3 = 
0.05 (left panels); 0.10 (right panels). 
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Figure 2: The same as in Fig. ^ but for yU~ and /i"*" event rates integrated over the neutrino 
(and muon) energy in the interval E = (2.0 — 20.0) GeV. 
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Figure 3: The same as in Fig. ^ but for fx and /i"*" event rates integrated over the neutrino 
(and muon) energy in the interval E = (2.0 — 20.0) GeV and |Am|^| = 3 x 10~^ eV^. 
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Figure 4: The same as in Fig. ^ but for /i and /i"*" event rates integrated over the neutrino 
(and muon) energy in the interval E = (5.0 — 20.0) GeV and for lAmg^^l = 3 x 10~^ eV^. 
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Figure 5: The charge asymmetry A^-f^+ of the multi-GeV /i~ and /i^ event rates, in- 
tegrated over the neutrino (and muon) energy in the interval E = (2.0 — 10.0) GeV 
and over the Nadir angle in the interval corresponding to 0.30 < cos^„ < 0.84 {man- 
tle bin), as function i) of sin^2^i3 for |Am|;^| = 2 x 10~^ eV^ (left panels), and ii) 
of |Am3i| for sin^ 26113 = 0.10 (right panels): ^^-^+ (solid lines), (AJ^-^+)nh (dashed 
lines) and (Aj^'l +)ih (dotted lines). The results shown are obtained for sin^^23 = 
0.36 (upper panels); 0.50 (middle panels); 0.64 (lower panels). 
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Figure 6: The same as in Fig. but for /x and /i^ event rates integrated over the Nadir 
angle in the interval corresponding to 0.84 < cos^„ < 1.00 {core bin). 
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Figure 7: The same as in Fig. but for | Ar 

integrated over the neutrino (and niuon) energy in the interval E = (2.0 — 20.0) GeV. 



25 



0.45 



0.4 - 

0.35 3;, 

0.3 - 

0.25 - 



0.2 - 



0.45 



1— r 



1— r 



J L 



J L 



0.4 - 
0.35 znlll 



0.3 - 

0.25 - 

0.2 - 



0.45 



TT" 



TT" 



J L 



J L 



0.4 h 
0.35 
0.3 



T" 



0.25 - 
0.2 - 



0.15 



TT" 



TT" 



"\ 



J L 



J L 



1— r 



J L 



TT 



J L 



TT 



J L 



1— r 



J L 



TT 



J L 



TT 



J L 



T- r 



J L 



TT 



J L 



TT 



J L 












0.03 0.06 0.09 0.12 0.15 1 1.5 2 2.5 3 3.5 4 



sin 29 



13 



IAm^3^l [10"^ eV^] 



Figure 8: The same as in Fig. but for | A?7i|^| = 3 x 10^'^ eV^ and /i~ and /i"*" event rates, 
integrated over the neutrino (and muon) energy in the interval E = (2.0 — 20.0) GeV and 
over the Nadir angle in the interval corresponding to 0.84 < cos6'„ < 1.00 {core bin). 
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Figure 9: The same as in Fig. but for /i and /i"*" event rates integrated over the neutrino 
(and muon) energy in the interval E = (5.0 — 20.0) GeV, and for |A?7i3^| = 3 x 10"'^ eV^. 
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